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G^ ■ Abstract 

O ■ A central object in the interpretation of quantum mechanics of closed 

'^ ' systems with decoherent histories is the decoherence matrix. But only for 

$-H , a very small number of models one is able to give explicit expressions for 

Olj' its elements. So numerical methods are required. Unfortunately the di- 

mensions of this matrices are usually very high, which makes also a direct 
numerical calculation impossible. A solution of this problem would be 
^\^ ' given by a method which only calculates the dominant matrix elements. 

^ . This includes to make a decision about the dominance of an element before 

it will be calculated. In this paper I will develop an algorithm that com- 
bines the numerical calculation of the elements of the decoherence matrix 
with a permanent estimation, so that finally the dominant elements will 
be calculated only. As an example I apply this procedure to the Caldeira- 
Leggett-modell. 
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1 Introduction 

The Copenhagen interpretation of quantum mechanics rests on the artificial 
division of reahty in quantum systems and measuring devices. The resuh of 
a measurement is produced by a collapsing wave function, which is induced 
by interaction with a measuring device. This is inadequate for a quantum 
cosmology, where the universe is to be considered as a quantum system. By 
definition, there is nothing else outside the universe dealing as a measuring 
device. An interpretation of a quantum theory of closed systems is needed. 
Such an interpretation has been developed by Griffiths [Q, Omns (2|, Gell- 
Mann, Hartle |^, Joos, Zeh |j] and others. It makes no reference to external 
observers, classical apparatus or the collapse of the wave function. 

The main feature of this formulation is to consider a time sequence of possible 
outcomes of measuring processes. They are called histories. Fine-grained or 
elementary histories are the basic set of histories, from which all other coarse- 
grained histories can be obtained by unification in the sense of set theory. The 
way, the coarse-graining is done depends on the physical properties of the system 
one is interested in. 

To connect the coarse-grained histories with observables, one has to assign 
each of these coarse-grained histories a probability. The problem is that some 
of these histories may have a quantum interference, which violates the sum 
rules of the probability theory. These interference's are expressed by means 
of a decoherence matrix. If all the non-diagonal elements of this matrix are 
vanishing, the histories will be decoherent. Only in this case one is able to 
assign a probability to each history. 

To find the coarse grained sets of decoherent histories is a central problem 
of this quantum theory of closed systems. One can define a decoherence matrix. 
Its elements describe the quantum interference between two histories. 

Given the histories, the decoherence matrix can be calculated easily in prin- 
ciple. The usual type of calculations is simple, but their number is very large, 
because the number of possible histories that is related to the numbers of ele- 
ments of the decoherence matrix is usually very huge. This makes it impossible, 
to calculate the decoherence matrix even numerically. For example: Consider a 
particle in one- dimensional space divided in m intervals. A coarse-grained his- 
tory should correspond to the situation, in which the particle has been detected 
in one of these intervals at the times (ti, ^2, ■ • ■ , tn)- The number of possible 
histories would be to" and the number of elements of the decoherence matrix 
would be ■m^^"-^^\ 

Not all is lost, we are not interested in the decoherence matrix for any 
coarse-graining. We are usually interested in such coarse-grainings, that provide 
decoherent histories. So if we combine the process of constructing such coarse- 
grainings together with the calculation of the decoherence matrix, we might be 
more successful. In this case for instance we already know, that the non-diagonal 
elements are vanishing. 



So we can try to calculate the decoherence matrix in two steps. In the first 
step we have to find qualitative features of the decoherence process to get the 
orders of magnitude for a coarse-graining that guaranties the dominance of the 
diagonal elements of the decoherence matrix. In a second step we can use these 
data for a numerical calculation. 

In this paper we want to show how to do these two steps, especially how to 
design the numerical calculation for the Caldeira - Leggett model, which gives 
us some deeper insight in the mechanism of decoherence. 

2 Qualitative features of the Caldeira-Leggett 
model 

This model has been developed by Caldeira and Leggett 0], to explain the 
Brownian motion. It consists of one distinguished harmonic oscillator in an 
interaction with a huge number of other harmonic oscillators, which we will call 
the environment. The action of the distinguish oscillator is 

Sa[x\^— / dtii' -u^x^) (2.1) 

2 Jo 

and the action of the oscillators of the environment are given by 






(2.2) 



The meaning of the variables is obvious. The interaction between the environ- 
ment and the distinguished oscillator is given by 



y^ / dtCkRkX 
,.. JO 



Si[x,R] = -2_^ dtCkRkX. (2.3) 

The Cfc's are coupling constants. At an initial time ti the density matrix of the 
whole system should be 

p{xi,Ri,yi, Qi;ti) = pA{xi,yi;ti) pe{Ri, Qi) (2.4) 

and we assume that the density matrix of the environment is in thermal equi- 
librium. 

Pe{Ro,Qo) = l[cxp[-j-^^jj^^[{Rl + Ql)cosh{u;k/kT)^2RkQk] 

2t: smh{ujk / kT) 

The space of all possible results of a measuring process is the phase space 
(x, i?;p, P) with positions x and momenta p of the distinguished oscillator and 



R, P are the coordinates and momenta of the environment. An elementary 
history corresponds to the classical situation, that at the times (ii, t2, ..., tn) 
the system can be found at {{xi,Ri;pi,Pi), (2:2, i?2;P2, ^2), •■•, {xn,Rn;Pn,Pn))- 
We are only interested in the position of the distinguished oscillator, neither in 
its momenta nor in the phase space coordinates of the environment. In this sense 
we collect all elementary histories with equal position coordinates 0:1,0:2, ■■■,Xn 
together in one family. The most general description of the above defined fam- 
ilies of histories is given by the reduced density matrix 

p{x, y;t)= f dRdQ6{R - Q)p{x, R, y, Q; t) (2.6) 

Now we want to investigate the time development of the reduced density matrix. 
Following the papers by Caldeira and Leggett Q,and Halliwell and Dowker 
[|6| we can represent the evolution of this density matrix with the help of a 
propagator J 



p(a^/>y/;*)= / da;odyoJ(a;/,y/;t|xo,yo;0) p(xo,yo;0). (2.7) 

By replacing the huge number of oscillators of the environment by a continuum 
with the distribution 

, ^„2^ ^ f 4Mm7W^/7r ,uj<Q. ,„ „, 

Pd{-)C\u)^^ \oj>n (2.8) 

(the constant 7 represents an effective coupling) the authors could give an ex- 
plicit formula for the propagator in the high temperature limit: 

J{xf,yf- t\xo, yo; 0) = F\t) exp(z5,, - $) (2.9) 

where the functions in the exponent are as follows: 

S,i = K,{t) {x) - yj) + Kit) {xl + j/2) 

- L{t)ixo+yo){xf-yf)-N{t)ixf+yf){xo-yo) (2.10) 

and 

<^ = A{t) {xf - yff + B (t) {xf - yf) {xo ~ y,) + C (t) {xo ~ yof (2.11) 

The explicit expressions for the time dependent functions Ki, K2, L, N and 
A, B, C will be given in the appendix, they consist mainly of the exponential 
function exp(— 7t) and trigonometric functions in t. This propagator has the 
remarkable property, that its time dependence includes only the difference be- 
tween the initial and final time. Originally one should have expected, that 
the propagator is not time-like translational symmetric, caused by a time de- 
pendence of the environment. But by replacing the discrete oscillators by a 



continuum, the environment became infinitely large and so the influence of the 
distinguished oscillator on it can be neglected. 

With the explicit formula (2.1Cl|) we have all information about the evolution 



of the density matrix. But the structure of these formulae is very complex 
and the high number of parameters makes it difficult to analyze the evolution 
in detail. So I decided, to study the decoherence matrix numerically. The 
result was, that for a sufficiently large final time we got asymptotically always 
the same density matrix. This reflects the fact, that the oscillator comes to a 
thermal equilibrium with the environment. Now we want to confirm this result 
analytically. 

For the moment we assume, the initial state of the distinguished oscillator 
at the time ii can be described by the Gaussian density matrix 

PA(a;i,2/i;i) = 1/ exp( ). (2.12) 

with arbitrary parameters a^ , Cy and a^y ■ This density matrix is normalized 
by the function V. The function pA represents a Gaussian wave packet in the 
"squared" configuration space of the distinguished oscillator. The variables ^o 
and 770 denote complex numbers, which give an information about the position 



and "momentum" of this packet. Using the formula (2.7) it is easy to get 
the density matrix at a later time. One has simply to perform a Gaussian 
integration. The result is again a Gaussian concerning the new position variables 
Xf and Hf. The coefficients are rational functions of the functions A, B,C, . . . ,N 
and thus rational functions of exponential and trigonometric functions. So one 
should assume, that these coefficients can be expanded in a power series with 
respect to exp(— 7*) for large t. These calculations are very tedious, but can be 
carried out by means of computers. The result is: 

+ exp(-7t) (a(i) Xf + b{t) yj) 

+ 0(exp(-27f))j (2.13) 

The functions a{t) and b{t) describe a harmonic oscillation with the frequency 
Lo. The evolution of the density matrix in the squared configuration space can 
be splitted in three phases: In the first phase, which has a characteristic period 
of r = "yt/2, the shape of the wave packet relaxes to a certain form, that 
does not depend on the initial parameters of the packet at ti. In the squared 
configuration space this form represents a two dimensional Gaussian bell which 
is located over an ellipsis iS. It corresponds to a certain " equi-density" line. 



This can be seen from our formula (2.13) by the fact, that all terms with the 



coefficient exp(— 7^) become of smaller order compared to others. In this sense 
they can be neglected and thus the coefficients of the terms quadratic in Xf and 



yf become independent of time and also independent of the initial parameters 
of the packet. From formula (2.13|) we can see also, that the main axes of the 



corresponding ellipsis S have a length of 

dB^ = "' 'SJ^'^ (2.14) 






and 



dg = ^— . (2.15) 

They are oriented at an angle of 45 degrees against the coordinate axes. 

While the shape of the packet becomes fixed after the first phase, its motion, 
corresponding to a damped oscillation of its position in squared configuration 
space is still present. This can be considered as to be the second phase. Its 
characteristic time is also r. After that time the system rests on the origin of 
our coordinate system in the squared configuration space. This can be con- 
sidered as the beginning of the third phase. In that phase the distinguished 
oscillator is in thermal equilibrium with the environment. If we would remove 
the coupling, e.g. set 7 = 0, then this oscillator will be on equal footing with the 
oscillators of the environment. The thermal equilibrium can then be recognized 



from our asymptotic formula (2.13) by the fact, that it is just the high temper- 



ature expansion of one of the factors in (2.5), but of course with the specific 
parameters of the distinguished oscillator. For 7^0 the whole calculation can 
be understood as the quantum analog of the damped harmonic oscillator. 

By considering, that all initial functions for the density matrix can be com- 
posed by a linear superposition of Gaussian waves (This system is even over 
complete), we see, that the asymptotic behavior, described above, is very gen- 
eral. 

Now we want to use these results, to construct a physically interesting set 
of decoherent histories and to make statements about the corresponding de- 
coherence matrix. At first we have to enlarge our coarse-graining. From the 
Heisenberg uncertainty relation we know that it makes no sense, to measure ex- 
act positions of the distinguished oscillator at the times ii, ^2, ■•■, in- We have to 
ask questions with a higher degree of coarse-graining. So we ask only, if the os- 
cillator goes through definite intervals at these given times. Our coarse-graining 
is now, that we collect in one family all histories, in which the position of the 
distinguished oscillator goes through a sequence of intervals A^^, A|^, ..., Af^ at 
the times ti,t2, ...,tn. In our notation for the intervals, the upper index means, 
that we may use different sets A'^' of intervals for different times t^ and the lower 
index denotes which interval of a set we selected. Sometimes we use also the 
notation 

A,„,,....,„=A,\,A2^,...,Ar„. (2.16) 

One family is now understood as one of the new coarse-grained histories. If 
the sets of intervals for one particular time are disjunct and their union covers 



the whole space of positions, the corresponding histories will be complete and 
alternative. 

The decoherence matrix corresponding to the coarse-graining defined above 
is given by 

^(A»i,...,^„A,,,...,,J = 

dxndynS{xn - y„) 

A" =A" 

^ '. Jn 

dXn-idyn-lJ{Xn, Vn, tn\Xn-l, Vn-i; Ui-i) 



dXn^2dyn~2J{Xn-l,yn-i;tn-l\Xn-2,yn-2;tn-2) 

-2 Jn-2 



A2 .A^ 

^2 ■ J2 



dx2dy2J{x3,y3;t3\x2,y2;t2) 
dxidyiJ{x2,y2; t2\xi,yi;ti)pAixi,yi;ti) 



Ai .Ai 



This describes a sequence of evolution's between two subsequent moments 
of time and following projections. The symbol Aj^ i„bi,...,i„ denotes the set 
of intervals A^^, ..., A"^, A^^, ..., A"^. We introduce two additional terms: The 
first is the incomplete decoherent matrix at time ife. Its elements are defined by 
the expression 

P-e(A,,,...,,,_,|j,,...j^_Jxfc,yfc,tfc) = 

dxk-idyk-iJ{xk,yk;tk\xk-i,yk-i;tk-i) 






(2.18) 



dxk-2dyk-2J{xk-i,yk-i;tk-i\xk-2,yk-2]tk-2) 

*fc-2 Jfc-2 

dx2dy2J{x3,y3;t3\x2,y2;t2) 

2 .A^ 
*2 ■ 32 

dxidyiJ{x2,y2;t2\xi,yi;ti)pA{xi,yi;ti) 

Ai ,Ai 



A second term is the projected incomplete decoherent matrix at time tk- The 



elements are given by the expression 



P+e{\i,...,ik^^jl,...Jk\Xk,yk]tk) 
P-e(Aji,...,j^,_i,A. 

, otherwise 



P^e{'^ii,...,z,^,,^ju...jk-i \xk,yk, tk) , XkeA'y^ , y^eA^^ 



(2.19) 
The physically interesting question is, for which sets of intervals are the cor- 
responding histories decoherent? In all what follow, we assume the timelike 
distance between two "measurements" is larger than the characteristic time r. 
The decoherence of the histories is equivalent to vanishing non-diagonal ele- 
ments of the decoherence matrix. Any of these elements corresponds to a pair 
of histories which have different intervals at least at one moment of time. Let 
this time to be tk- The incomplete decoherence matrix p-^(tk) can be written 
as 

P-e{A^U...,^k^^ln,■■■,3k-Mk,yk■,tk) = 

dxk-idyk-iJ{xk,yk;tk\xk-i,yk-i;tk^i) 

^P+<i{^^l,...,^k-l\3l,...,jk-Mk-l,yk-l]tk-l) (2.20) 

This is an evolution of p^^{tk^i) from tk-i to tk- If this time interval is larger 
than the characteristic time r, we can assume that p_^{tk) is equal to the 
asymptotic matrix pg. The subsequent projection of the incomplete decoher- 
ence matrix on the intervals A*^ and A*^ vanishes approximately, when these 
intervals do not catch the ellipses. This is guarantied when the distance of the 
centers of the intervals will be larger than the width ds and because these two 
intervals must be different and disjunct this is given, when their own width is 
larger than that of ds- 

If the incomplete decoherence matrix p^^[tk) vanishes, then the chain of 



evolution's and projections in formula (2.17) is cut off and the corresponding 



matrix element vanishes also. Because of the criterion, that the intervals for 
projections should be larger than ds does not depend on any other special 
property of the histories we considered, it guaranties decoherence of all histories 
with that property. So we will call ds the decoherence width from now on. A 
remarkable fact is, that this width does not depend on 7, measuring the strength 
of the coupling. But the timescale for the decoherence does. 

In what follows we will always assume, that the widths of the projections 
are larger then the decoherence width. So our histories will decohere and we 
ask now for their probability. For calculating the probabilities the projection 
intervals Aj"' and A*^ have to be identical. We work in the configuration space 
of the distinguished oscillator. 

Again we start with the incomplete decoherence matrix p^^{tk), which comes 
from a projection at tk-i and has the asymptotic form of a wave packet ps- 



Independent on which interval we project at tk, the asymptotic form, at tk+i 
will be the same. Only the prefactors will become different. Are the intervals 
for projection at tk outside the wave packet of p-eitk), the prefactor at tk+i 
will be zero. Otherwise it will always have approximately the same order of 
magnitude if it goes through the packet. This means that we get approximately 
the same probability for all histories, whose fc-th interval is completely within 
a region, which is the projection of the main axis with dg of the ellipsis S to 
the coordinate axis x. The width of that projection is 



'^ " kT 



dB = " .^^ (2.21) 



If we apply the same argument for all other times, we see that all histories 
having intervals only inside a strip of width ds in the configuration space of the 
distinguished oscillator, have the same probability. Is the with ds larger than 
the asymptotic decoherence width ds, then the oscillator shows a Brownian 
motion within an interval of ds ■ 

By analyzing the decoherence width ds and the width of the Brownian mo- 
tion we find the familiar result, that with increasing temperature T the decoher- 
ence length becomes smaller, while the Brownian motion becomes more intense. 
A stronger coupling makes the Brownian motion weaker, but has no infiuence 
on the decoherence length. A bigger mass makes the decoherence length and 
the Brownian motion smaller. 

This calculation was for the Caldeira-Leggett model. But the fact, that in 
the process of evolution the distinguished oscillator comes into thermal equi- 
librium with its environment in general, makes it very likely that in all similar 
calculations we will find a specific asymptotic behavior of the observed subsys- 
tem. This might give the possibility, to set up a quite similar calculation to get 
the decoherent histories for other more general systems. 

It is also interesting to investigate the evolution of the density matrix for 



short times. For our initial density matrix ( 2.12 ) the width for decoherence is: 



do' -i^ + ^- -^)/4 (2.22) 

Ox Uy (J XV 



An expansion of the density matrix for small t gives the following result: 
Ps{xf,yf;t) = y'exp 



i^f-yff{^ + i'2kTM-^)^t + 



(- 



* 2Malala,y J 

f + yf\2({(Tx(^xy + '^y<^xy + <^x(^y){Max<Jy + '2.i{ax ~ Cry)t)-' 



Ma^a^axy 



+ CxXf + Cyjjf + Co 



(2.23) 
The variables C'x, Cy and Co represent some time dependent function in which 
we are not further interested. We focus our attention on the real part of the 
coefficient of {xf — y/)^. It describes the evolution of the decoherence width D. 
If we replace the temperature by means of the asymptotic decoherence width 
ds, see ( ^.15 ), we get 



1 
15^ 



1 



-4( 



1 

4 



1 
-4 



)it 



(2.24) 



We see, is the initial decoherence width smaller or larger then the asymptotic 
one, the width will be immediately increased or decreased respectively. So 
we can say that the relaxation starts immediately. Suppose the initial width 
do is much larger then the asymptotic width dg. How long does it take, till 
the decoherence width of the oscillator is smaller than a given value of d. By 
assum ing t hat this width is also much larger than the asymptotic value, we get 
from (12^2^) 

Depending on the ratio d^/d t his ti me might be much shorter, than the relax- 
ation time T. From formula ( 2.23 ) we see, that the time dependent part of 
the width of the Brownian motion is imaginary. So the range of the Brownian 
motion is not influenced by the time dependence in first order. This provides 
another scenario than described above. If we consider coarse-grainings with in- 
tervals much larger than the asymptotic decoherence width, then in a first phase 
of the scenario the desired decoherence happens in a very short time compared 
to the relaxation time. After that the damped oscillation and Brownian motion, 
known from the asymptotic scenario, takes place. 

If we choose a coarse-graining with intervals in the order of magnitude of 
ds and the time between measurements in the order of r, the diagonal and 
"near-to-diagonal" elements of the decoherence matrix should be dominant. So 
we can neglect the others and the number of elements we have to calculate is 
numerically manageable. 
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3 How to calculate the decoherence matrix nu- 
merically 

In the case of Hamiltonian quantum mech anics the decoherence matrix can be 
calculated by a sequence of integration's ( 2.17 ). These integration's might be 
carried out numerically. But the number of matrix elements is very large in 
practical cases. So, to get reliable numerical results one has to find out the 
dominant matrix elements, without spending to much time for non dominant 
terms. Therefore we want to elaborate some inequalities, to give upper bounds 
to the matrix elements. 
From ( ^.20 ) we get 



D{A,, 



,in\jlr---Jn 



/ dXn-ldyn-1 

UXjiJ \Xji^ ^m ^n l-^n— 1 7 Vn—l 1 ^n— 1 / ^ 



(3.26) 



By Cauchy's inequality with respect to the integration's over dxn-idyn-i and 
with the following definitions: 



a(a^jaj:-;,a;:::^)^ = 






dxn-idyn-1 



A" 



ClXjiJ {^Xfi , Xn , til I'^n — l j yn— 1 1 ''n— 1 ) 



and 



[,..,h\3k,..,ji = / dxkdyk\p-eiA,^,..,r,\j„...,jJxk,yk]tk)\ 

J A'' .A*" 



(3.27) 



(3.28) 



we obtain from (3.26) the inequality 

|i^(A,,,...,„|,,,...,,J| < A(A^IA-:\,A-:^)J^„_,,..,,,|,„_,,..,,, (3.29) 
Now, we want to find an upper bound for Pik,..,ii\jk,--,ji ■ With the help of formula 
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(^.2C) we can rewrite equation (B.28) and after rearranging the integrals we get 









»fc-ll3fc- 



dxkdykJ{xk,yk;tk\xk-i,yk^i;tk-i)J{xk,yk;tk\x'k_i,yk-i;tk-iy 



Again by using Cauchy's inequality and with the definition 
i3^AtA^|At^A^)=/ dxk-.dyk-J dx',,,dy',^. 



(3.30) 



»fc-ll3fc-l 



dxkdykJixk,yk;tk\xk-i,yk-i;tk-i)J{xk,yk;tk\x'^_i,y'k-i;tk-iy 



(3.31) 



we can write 



P, 



»fc,--,2i|jfc.-- J 



. < 5(AtA^|At^At^)i^,_,,..,,|,,_,,..,,. (3.32) 



If we know the incomplete density matrix at a t ime slice tk, w e can calculate 
Pik,..,ii\jk,..,ji ^-nd an iterated use of formula ( ^.32 ) together with (|3.29 ) gives the 
possibility to give an upper bound for the elements of the decoherence matrix at 
any later time. This is not yet very useful for a numerical calculation, because we 
can not test the upper bound for all possible path, caused by their large number. 
So we have to ask: suppose we know an incomplete density matrix at a time 
tk, what is the upper bound for all decoherence matrix elements independent of 
subsequent projections? To answer this question we introduce the quantities 



A„,„_i(A^^-;,AJ-; 

^n.n— 1 



max A" 



max^„ i^.-i ^,^-1 



Bk,k-i{Al_\,A';Jj =maxA. ^A'- 



B. 



k,k-l 



max^fc ^k ^^-^ A^^^ 



With these maxima we construct the quantity 



,(fc) 



(3.33) 



^il',..Mik,..,jl = '^n,n~lBn-i^n-2,--;Bk+l,k{A-^,AjJP,^^,,^,^y^^,,j^ . (3.34) 
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The main point is, that the elements of the decoherence matrix satisfy the 
inequahty 

D(A,, , I,- , ) < K^'"'^ . ,. . . (3.35) 

It says that aU matrix elements containing the intervals Aj,,^ j^Uj. ^^ are smaller 
than the corresponding value of K^^\ independent of subsequent projections. 
Now we can use a computer in the following way: We start with the initial 
density matrix pA{xi,yi) and calculate the quantities 

Pn\n =" / dxidy2 \pA{xi,yi)\'^ (3.36) 

and K {., . We denote with M the set of all K , ., 's and consider the maximum 
«iui nbi 

of its elements to be the norm of A4 

M = ^nlniKlllJ (3.37) 

IIA^II = max(i^(;)^.J (3.38) 

Now we take the largest element oi M, the attached intervals should be A} Aj 
and propagate the corresponding incomplete density matrix. This allows us to 
calculate the quantities 

PlMn.h= i dx2dy2\pU^i,\.j,\{x2,y2-M)\' (3.39) 



,A^ 



and K^ ^ , . j . We enlarge our set M. by the quantities K'^'^^ and reduce it by 

the element Kj L . These steps will be repeated: We select the largest element 
of A4, propagate the corresponding incomplete density matrix and replace this 
largest element of A4 by the new obtained K^s. So for instance, we could have 
after two steps one of the iiT'^^'s as the largest clement or one of the iiT^^^'s. 

By repetition of these steps one of our propagation's reach the final timeslice 
tn and we can calculate the corresponding element of the decoherence matrix. It 
is not guarantied that this element is the largest, but we know that all elements 
we have not yet calculated are smaller than the norm of Ai . By going on with 
our procedure we get more and more matrix elements while the norm of A4 is 
decreasing. If this norm becomes smaller than some desired order of magnitude 
of the already calculated matrix elements, we can neglect the not jet calculated 
elements and the problem is solved. The only bad thing that can happen is, 
that the histories are not decoherent enough. So the program will produce more 
and more matrix elements, but the norm of yVf will drop down very slowly. If 
this happens one can find from the calculated data, which histories are not 
decoherent and one can combine some of the intervals to increase the coarse- 
graining. 

13 



To demonstrate the efficiency of the algorithm we consider the Caldeira- 
Leggett modeL The decoherence matrix should correspond to the situation 
where we have 5 intervals at 6 time slices. The parameters of the model and 
the intervals are chosen in such a way, that the widths of the intervals are 
one asymptotic decoherence width ds and the distance of the time slices is 
half of the relaxation time. With these parameters the decoherence matrix has 
5 X 5^ X 5^ « 4.8 X 10^ elements and for each element we have to calculate 
5 propagation's of the incomplete density matrices. The time to calculate the 
decoherence matrix is proportional to the number of projections 

P « 2.4 X 10^ (3.40) 

Let us assume, that we are not interested in matrix elements that are smaller 
than 10^^ times the largest element. 

Furthermore, for the numerical calculation we arrange the parameters of the 
model such, that the mass of the oscillator corresponds to that of a proton, 
the asymptotic decoherence width should be equal to one atomic length scale 
L = 10~^'^'m and the asymptotic width of the Brownian motion will be 3L. 

I made the additional assumption, that the oscillator is in a pure state at 
the initial time and that the wave function is a Gaussian function with a width 
of 2L, located at the origin of the coordinate system. 

The program started to calculate 8478 elements of the decoherence matrix. 
Whenever the program found that the corresponding matrix element will be- 
come smaller than one tenth of the largest element, it terminated that calcu- 
lation, so that finally only 1536 relevant elements had been determined. The 
number of propagation's, which had to be performed was 6025. The ratio be- 
tween the number of executed propagation's and the total number P is a measure 
for the efficiency of the algorithm and has a value of 2 x 10~^. This numerical 
calculation, should simply be considered as a demonstration. 

4 Discussion 

The purpose of this paper has been to provide a guide line for calculating the 
elements of the decoherence matrix. I developed an algorithm which combines 
the numerical calculation of the elements of the decoherence matrix with a per- 
manent estimation, so that finally the dominant elements will be calculated 
only. The algorithm can not be used straight forward to calculate other mod- 
els. Because other models require different inequalities to estimate the matrix 
elements. To find these inequalities may require many analytical calculations 
- better estimations provide a more efficient algorithm because nondominant 
matrix elements will be ruled out earlier. This has two advantages: the number 
of steps to be performed is smaller and the amount of data to be saved tem- 
porarily. But despite these model-dependent inequalities the main features of 
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the algorithm remain the same, so that paper may help to study more complex 
models. 



A Appendix 

The explicit formulae for the time-dependent Coefficients in the exponent of the 
propagator J ( 2.10| ) and ( 2.11 ) are as follows: 



kTM , , 2 7^ + t^^ - 7^ cos(2wt) - 7^ sin(2wT) - e-'^^^u? 
esc (wr) 



2 

h = kTMui csc{lu t 
kTM 



■ csc{ijJtY 









72+C^2 






2 


LO COs{ijJt) (1 — 


g2 7r^ +7 sin(a;T 


)(l + e27 











eT^ (72 + ^2) 






' 


-1'- 


u? + e2^- 


u? +72cos(2cjt) 


— 7a;sin(2 


lot) 



72 + a;2 



and 



Af7 Mbj , , 

A/7 Mw , , 
K2 = +^ + ^-cot(^r) 

. Mu _^^ 
L = -e ' 



N = 



2 sin (ojt) 



2 sin (ojt) 
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